Rapidity and kr dependence of HBT correlations in 
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Abstract. Two-particle correlations of identical charged pion pairs from Au+Au 
collisions at y/s^ = 200 GeV were measured by the PHOBOS experiment at RHIC. 
Data for the most central (0-15%) events were analyzed with Bertsch-Pratt (BP) 
and Yano-Koonin-Podgoretskii (YKP) parameterizations using pairs with rapidities of 
0.4 < y < 1.3 and transverse momenta 0.1 < kx < 1.4 GeV/c. The Bertsch-Pratt radii 
decrease as a function of pair transverse momentum. The pair rapidity Y^^r roughly 
scales with the source rapidity Yykp, indicating strong dynamical correlations. 

Identical-particle correlation measurements (Hanbury-Brown and Twiss, HBT) yield 
valuable information on the size, shape, duration, and spatiotemporal evolution of the emission 
source in heavy ion collisions. Experimentally, the correlation function C(q) is defined as 

^^^^ " P(Pl)P(P2) ^'^ 

where -P(pi,P2) is the probability of a pair being detected with relative four-momentum 
q = pi — p2, and -P(pi) and -P(p2) are the single particle probabilities. The numerator 
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is determined directly from data, while the denominator is constructed using the standard 
event-mixing technique. 

The data reported here were collected using the PHOBOS two-arm magnetic spectrometer 
during RHIC Run II (2001). Details of the setup have been previously described in pp. The 
spectrometer arms are each equipped with 16 layers of silicon sensors, providing charged 
particle reconstruction both outside and inside a 2 T magnetic field. The primary event trigger 
was provided by two sets of 16 scintillator paddle counters, which covered a pseudorapidity 
range 3 < |r/| < 4.5. Details of event selection and centrality determination can be found in 
[21 E]- The 0-15% most central events were used in this analysis, equivalent to (Npart) = 310 
as determined by a Glauber model. 

The details of the track reconstruction algorithm can be found in [1]. Events with a 
reconstructed primary vertex position between -12 cm < z^tx < 10 cm along the beam direction 
were selected in order to optimize vertex-finding precision, track reconstruction efficiency, 
and momentum resolution. Only particles which traversed the entire spectrometer were used 
in the analysis. A 3a cut on the distance of closest approach with respect to the primary 
vertex (dca^tx < 0.35 cm) was then applied. The final track selection was based on the x^ 
probability of a full track fit, taking into account multiple scattering and energy loss. The 
momentum resolution is Ap/p r^ 1% after all cuts. To identify pions, a cut three RMS 
deviations away from the expected mean value of the specific ionization {dE/dx) for pions was 
applied. Contamination from other particle species was studied using HIJING 1.35"^ and a 
GEANT 3.21 simulation of the full detector. The contamination from K K , pp, and p p pairs 
is less than 1%; non-identical pairs contribute less than 10% throughout the entire /ct range. 
To reject ghost pairs, only one shared hit in the weak-field region and two shared hits in the 
strong-field region were allowed per pair. A two-particle acceptance cut was applied to both 
data and background; the criterion for pair acceptance was defined by Ac/) -|- 2A0 > 0.05 rad, 
where A0 and /S.6 are the relative pair separation in azimuthal and polar angle, respectively. 
About 7.3 million Tr+vr^ and 5.5 million vr^vr^ pairs survive all cuts. 

Systematic errors were determined by changing two-particle acceptance cuts, cuts in 
azimuthal separation, random seeds used in mixed-event background generation, as well as 
varying the definition of "event class" to create background events from pairs within narrow 
and broad vertex ranges. 

Because the event-mixed background is the product of tracks from different events, it does 
not o priori include any multiparticle correlations. In order to study the HBT correlation, it 
is necessary to apply a weight to account for the Coulomb effect. The Coulomb correction can 
be expressed solely as a function of relative 4-momentum q, 

^^^^ FM Jdr\i;,i{r)\^S{r) ^' 

where S{r) is the relative separation of the particle pair, and ipc and ippi are the Coulomb and 
plane wave- functions, respectively. A closed-form approximation and numerical correction for 
this relation was derived in [£] for A = 1. For a variable A, 

.,,,.,.(l^^)±^ll±5p^ (3) 

This prescription is nearly equivalent to the corrections applied by the CERES, STAR, 
and PHENIX experiments [Zl El Eli our results showed no significant change using either 
correction method. The method is applied iteratively, successively fitting distributions of the 
correlation function C{q) and iteratively applying the fit value R to a. new S{f). Typically 2 
or 3 iterations are sufficient for convergence. 
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C(q) is typically fit to a Gaussian source in three dimensions, the so-called Bertsch-Pratt 
parameterization jlUj . 

C(q) = 1 + Xe'('ioRlWsRl+q!RJ+2qoqeRlt) (4) 

The correlation function was also fit to the YKP parameterization jllj . 

C7(q) = 1 + Xe-'^ilRl+yHn-i3ir?Rl+jHir^PqnfRl) ^5^ 



where P is the longitudinal velocity of the source and 7 = l/\/l — /3^, q± and gy the relative 
3-momentum difference projected in the transverse and longitudinal directions respectively, 
and Qr the relative difference in energy. In order to compare with lower energy, the data 
presented was fit in the longitudinal co-moving system (LCMS) frame. 
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Figure 1. Bertsch-Pratt radii as a function of pair transverse momentum fc^ for 
Au+Au at 200 GeV from PHOBOS, STAR 'H' and PHENIX ^. The boxes represent 
PHOBOS systematic error. 



In Fig. ^ the Bertsch-Pratt radii are presented as a function of pair transverse momentum 
kx for TT~TT~ pairs. For comparison, data from STAR !8| and PHENIX "W at Ws^j^ = 200 GeV 
are also shown. The PHOBOS data were analyzed in the LCMS frame within the rapidity 
range 0.4 < y < 1.3, while the other data are at mid-rapidity (—0.5 < y < 0.5). The three- 
dimensional correlation functions were fit to Eq. Q) using the log- likelihood method. Rg 
weakly varies as a function of kx, while Rq and Ri decrease rapidly with increasing kx- 

In Fig. 121 the extracted value of the source rapidity Yykp is plotted as a function of 
pair rapidity for tt+tt"'" pairs with 0.1 < kx < 1.4 GeV/c. The data from NA49 p^ at 
lower energy is also plotted; however, it should be noted the presented NA49 data covers only 
0.1 < A;^ < 0.2 GeV/c. The pair rapidity strongly scales with source rapidity, indicating the 
presence of strong position-momentum correlations. The solid line at Yykp = Y-k-k represents 
a class of models including, but not limited to, boost invariance. 
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Figure 2. Source rapidity (Yykp) as a function of pair rapidity {Y^.^) for PHOBOS 
(circles) and NA49 (squares) ^3. The line at Yykp = ^tttt is drawn to guide the eye. 
The boxes represent PHOBOS systematic error. 



In conclusion, we have extracted HBT parameters from Au+Au collisions at Ws^j^ = 
200 using two different parameterizations of the correlation function. The Bertsch-Pratt 
parameters show good agreement between three experiments with very different acceptances. 
Prom the YKP analysis, the pair rapidity scales strongly with the source rapidity, indicating 
a source with strong position-momentum correlations. 
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